How does the chromatin fiber deal with topological constraints? 
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In the nuclei of eukaryotic cells, DNA is packaged through several levels of compaction in an orderly re- 
trievable way that enables the correct regulation of gene expression. The functional dynamics of this assembly 
involves the unwinding of the so-called 30 nm chromatin fiber and accordingly imposes strong topological con- 
straints. We present a general method for computing both the twist and the writhe of any winding pattern. An 
explicit derivation is implemented for the chromatin fiber which provides the linking number of DNA in eukary- 
otic chromosomes. We show that there exists one and only one unwinding path which satisfies both topological 
and mechanical constraints that DNA has to deal with during condensation/decondensation processes. 

PACS numbers: 87.16.Sr 



In the cells of higher eukaryotes, e.g. animals or plants, 
meters of DNA are packaged by means of proteins into a nu- 
cleus of a few micrometer diameter, providing an extreme 
level of compaction. Coding sequences (genes) are there- 
fore dispersed in a mess of folded DNA and proteins (chro- 
matin) and should be retrieved at will in order to enable a cor- 
rect genetic expression and therefore the cell survival. This 
leads to the need of an orderly and dynamically retrievable 
structure, which is actually achieved by means of a chromatin 
partition into functional domains, each containing one or a 
group of genes. In each domain DNA is folded in a hierar- 
chical structure, including several winding levels. It is first 
wrapped around spools of proteins thus forming a "beads on a 
string" assembly, which is in turn folded into a 30 nm diameter 
fiber. This fiber is further organised into a three dimensional 
cross-linked network 1 1 ] . In this network, two neighboring 
nodes are connected by a chromatin fiber loop whose typical 
length is about 50000 base pairs (bp). In order to provide the 
transcription machinery with an access to specific genomic 
regions, the corresponding loop has to be selectively decon- 
densed, via a reversible unwinding process that elongates the 
fiber 1 2]. The dynamics of this process involves strong me- 
chanical and topological constraints, the former due to DNA 
elasticity |3], the latter due to the conservation of the linking 
number Lk of DNA |4] in a loop during the unwinding. Of 
course, topological constraints could be released by the inter- 
vention of topoisomerases, but it has been shown in vivo that 
chromatin decondensation could take place even when those 
enzymes were inhibited |5]. Moreover, the classic experiment 
on simian virus 40 (SV40) minichromosomes clearly demon- 
strates that the linking number is unaffected by the deconden- 
sation process in vitro |6]. 

In this Letter we address the issue of how Evolution has 
dealt with the extremely difficult problem of finding an effi- 
cient winding pattern fulfilling both mechanical and topologi- 
cal constraints at a time. To answer this question, we start by 
giving an analytical formula for the linking number of DNA 
in a generic bent fiber. We show that, despite Lk is known 
to be a non extensive quantity, it is yet possible to express it 
in terms of the mean linking number of the constitutive ele- 
ments of the "beads on a string" assembly. This allows us to 



set about an exhaustive numerical exploration of Lk for all the 
possible fiber conformations. By analyzing the results of this 
exploration, we are able to infer the existence and the unique- 
ness of a relevant winding/unwinding path satisfying all the 
constraints. We show furthermore that this engineering prob- 
lem is solved at the local level by the design of the constitutive 
elements of the fiber. 
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FIG. 1: Schematic representation of the DNA winding pattern along 
two neighbouring nucleosomes in the two-angle model, (a): view 
down the NCP axis a,-, (b): view down the linker direction t[l\ . The 
angle a,- is the dihedral angle (til\,ai,7i) and (3,_i the dihedral angle 
(ajll ,tili,3i) standing for the twist (modulo 2%) of the DNA linker. 
We also indicate the DNA radius r ~ 1.0 nm, NCP radius R ~ 5.3 
nm and NCP pitch P ~ 2.4 nm. 

In a typical chromatin fiber, the beads, called nucleosome 
core particles (NCP), are spaced at intervals of ~ 200 bp. 
Each NCP contains ~ 150 DNA bp, forming 1 \ turns of a left- 
handed superhelix 0], and two neighboring beads are con- 
nected by ~ 50 bp stretches of DNA, called linkers. The unit 
of a NCP and a linker is called a nucleosome. The number of 
DNA bp in a nucleosome is known as the repeat length. In 
order to describe the structural polymorphism of a chromatin 
fiber loop, we use the original two-angle model (Fig.^ intro- 
duced by Woodcock et al. 1 8]. In this model, where the linkers 
are assumed to be straight, the geometry of a fiber made up of 
N nucleosomes is fully characterized by two sets of angles 
Uj and P, (i = 1 . . .AO, specific of the NCPs and the linkers 
respectively Varying a, corresponds to changing the 

number of bp wrapped in the ft 1 NCP, whereas varying p, cor- 
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responds to twisting the i linker. 

In vivo, the fiber is most probably homogeneous within a 
given loop full , i.e. that the distribution of a, (resp. (3,) is 
peaked around a mean value a (resp. p). During the cell cycle, 
the biological activity is tuned by biophysical 1 12] and bio- 
chemical 1 13] parameters, through a variation of these mean 
values. In turn, a variation of a or p results in a change of the 
fiber internal structure and hence affects its degree of com- 
paction. This conformational change is necessarily accompa- 
nied by a change of the writhe Wr and the twist Tw of the DNA 
double helix. Consequently, the DNA linking number is a pri- 
ori bound to change according to the White-Fuller theorem 
Lk = Tw + Wr I14lll5ll . This poses the rather important ques- 
tion of whether and how the angles a and p can be changed, 
while maintaining a fixed DNA linking number. To answer to 
this question, we set about computing the DNA linking num- 
ber in a fiber as a function of the sets 0C; and P,-. This involves 
the computation of both the twist and the writhe of the DNA 
double helix in the fiber. 

Concerning the computation of the writhe, Fuller's first the- 
orem 1 16] states that, for a closed curve r(s) with s S [0,L], the 
writhe is related to the signed area A enclosed by the tangent 
indicatrix T(s), namely Wr = A/271— 1 (modulo 2). T(s) is 
the curve on the unit sphere traced out by the tangent vector 
7(s) to the curve r(s). Fuller's second theorem 11611 permits 
us to get rid of the congruence modulo 2 by computing Wr 
from the area S swept out by the unique shortest geodesic arc 
from an arbitrarily fixed point C on the unit sphere to the run- 
ning point T(s), namely Wr = S/2k. However, the calculation 
of the writhe of an open curve whose initial and final tangent 
vectors 7(0) and 7(L) do not coincide is rather subtle. We fol- 
low the procedure recently outlined by Maggs 11711 to obtain a 
consistent measure for the writhe by closing the open tangent 
indicatrix with a geodesic arc. This being rather complicated 
for a general fiber conformation, it would be helpful if one 
could calculate the writhe from the writhes of the individual 
nucleosome units. However, the writhe is known to be non 
extensive and, as shown by Starostin 1 18], the total writhe of 
a curve is given by the sum of the writhes of its parts plus 
the surface of the spherical polygon composed by the set of 
geodesies closing each part. In the case of a chromatin fiber, 
this non-extensivity may be overcome in the following way. 
As shown schematically in Figure|2] the DNA indicatrix T(s) 
can be naturally divided into N parts, namely the N nucleo- 
somes, in correspondence to % = 7(si), the direction of the ft 1 
linker. We denote by 7} the point on the unit sphere that cor- 
responds to fj. The total writhe of the curve, Wr, is then given 
by the following addition rule: 

Wr=^Wr i + ^f^. (1) 

!=1 Z% 

Wri is the writhe of the ft 1 nucleosome and is given by the sur- 
face S enclosed by the restriction of the tangent indicatrix be- 
tween 7f_i and 7} and the geodesic 7}_i7}. The non-extensive 
correction St Ti ...t n is the signed area of the spherical polygon 




FIG. 2: Schematic picture of the tangent indicatrix of DNA in a chro- 
matin fiber with the indication of the particular points used in the 
calculations: each solid-line arc corresponds to a nucleosome, with 
linker directions given by the connecting points 7). Points F,- indi- 
cate the local fiber axis directions. C is an arbirtarily fixed point. All 
dashed lines indicate geodesic arcs. For the sake of clarity, scales 
and lengths have been arbitrarily chosen. 



connecting all points 7] that correspond to the linker directions 
ti. It can be computed as the sum of the areas of the spheri- 
cal triangles 7b7/7} + i. It can also be expressed with respect to 
any point C on the unit sphere as 1 Sct^ i t, plus a boundary 
term Sct n t - Up to this term, the total writhe can be written as 
Wr = '£f =l Wr i \ c , where 

Wr i \ c = Wr i + ^^ (2) 

27C 

is the writhe of nucleosome i with respect to point C, cor- 
responding to the writhe of the nucleosome i alone plus the 
surface of the spherical triangle C7}_ 1 7} . We see that an ap- 
propriate redefinition of the writhe of a nucleosome permits 
us to express the total writhe as an extensive quantity. 

For the simplest case of a straight regular fiber all the nu- 
cleosomes are evenly dispersed on a regular helix, with the 
angles a, = a and p, = p Vz. In this case, though Wri is 
the same for all nucleosomes, Wri\c varies in general from 
nucleosome to nucleosome. Nonetheless there exists a spe- 
cial point C = F, defined by the director of the fiber axis, 
for which Wr^p is independent of i, so that we can define an 
effective writhe per nucleosome Wr(a, p) = Wr^p, such that 
Wr = NWr(a,$) UM- We reiterate that the effective writhe 
per nucleosome Wr(u, p) is not the bare writhe of one individ- 
ual nucleosome, the difference being the surface of the spher- 
ical triangle FTqT\ . 

For a generic bent fiber, a local fiber axis Fj can be defined 
for each i as the axis of the straight fiber characterized by con- 
stant angles a = a, and p = p,- and containing the ft 1 nucle- 
osome. As shown in Fig. this allows us to subdivide the 
surface of the spherical polygon TqT\...Tn into spherical tri- 
angles defined with respect to the points Fi, so that the total 
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FIG. 3: (a): solid contour lines (spaced at intervals of 0.05) represent the linking number Lk(a, (3) computed for ft ranging from 0° to 180°, (3 
from 0° to 360° and h = 192 bp. Color scale gives the fiber compaction defined as the number of nucleosomes per unit volume C = 1/ndR 2 , 
in units of NCP/ (30 nm)l Black regions correspond to sterically forbidden regions. Points 1 to 4 are particular reference states: structure 1 
is the most compact fiber structure allowed, paths 1-2, 1-3 and 1-4 correspond respectively to decondensation at constant (3, constant ft and 
constant Lk(&, (3). (b): fiber structures corresponding to points 1 to 4 in (a), top view and frontal view drawn for 24 nucleosomes. For clarity 
the first and the last NCP orientations are outlined by a yellow and a dashed violet line. The twist change between any two fiber structures is 
the angle swept out by the dashed violet line (namely, —0.70 between 1 and 2, —1.10 between 1 et 3, between 1 and 4). 
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up to the boundary terms 



writhe finally reads 

N 

with WV(a,,p,) = Wn 

Scf n f, +Sf n t n t f v 

As concerns the twist of the DNA double helix in the 
fiber, its extensivity enables a more straightforward calcula- 
tion. First of all note that, if the repeat length n, is allowed 
to change from nucleosome to nucleosome, the i™ linker will 
be relaxed (untwisted) only for a particular equilibrium value 
P^oc,-,/!,-) of (3,. Then, let rw(a;,P;) be the twist of the 
ft 1 nucleosome, with n, as a repeat length. At fixed a ; -, the 
variation A7V((X;, p,) with respect to p, is equal to Ap,/27t, 
whereas, at fixed p,, the variation with a, is equal to the varia- 
tion of the tortuosity xRAai/ln, where X = 2%P/(A% 2 R 2 + P 2 ) 
is the torsion of the double helix axis in the NCP i2^l . In 
order to have an expression for the total twist we just need 
to know, therefore, its value on a particular reference state. 
We take as a reference a nucleosome with 146 bp in the 
NCP 13] and a relaxed linker. In this state, a, = a ~ 60°, 
and p,- = P° = p"'(a ,n ; ) = (n; - 146)/% The twist is then 
Tw(a°, P°) = 146//; + (m - 146) /h where h (resp. h) is the 
helical repeat of the double helix free in solution (resp. over- 
twisted in the NCP) 12311 . The total twist of the DNA double 



helix in the fiber is therefore 
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For the DNA double helix in a straight and relaxed state, 
the writhe is zero and the total twist simply amounts to the 
last term in Eq. g), hence Lk° = Tw° = "i /ho =Nn/h , 
where h is the mean repeat length of the fiber. Therefore, 
defining the local twist with respect to the relaxed DNA as 
TV^Pi) - (T/?(a i -a ) + p ! -p?)/27t+146(l// ; -l//io), 
we obtain from Eq. and Eq. @ the following expression 
for the DNA excess linking number of a generic fiber: 
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is the local DNA 



where L/t(a,,P,) = Wr(a i ,p I ) + 7w(oc;,P; 
linking number per nucleosome. 

It is interesting to note that the very last term of Eq. 0, 
Yli=i ScFjF i+1 /2tc, is the area S swept out by the geodesic arc 
from C to the point Fi as i increases. Therefore, this term 
is nothing but the writhe of the fiber axis, Wr T . The excess 
linking number Lk — Lk° is the contribution of the winding 
pattern of the double helix to the total linking number Lk. It 
should therefore be interpreted as the linking number of the 
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fiber, Lk' f . The difference Lk T — Wr' 1 must in turn be identified 
with the twist of the fiber around its axis, Tw r , in order to 
satisfy the White-Fuller theorem at the level of the fiber Lk T = 
Wr' 1 + Tw T . This formal definition of Tw' f actually matches 
the intuitive definition of the fiber twist as the rotation angle of 
the fiber "top" with respect to the fiber "bottom", as shown in 
Fig. 13b) where we give the twist variations between different 
fibers. 

We now evaluate the relative contributions of the three 
terms in Eq. l|5}- The first one can be evaluated in a mean 
field approximation as Lfc(oc,-, P,) ~ NLk(a, p). We have 
therefore computed Lk(a, p) for a ranging from 0° (corre- 
sponding to "open" NCP wrapped with 1.5 turns of superhelix 
IU9I0 to 180° (corresponding to "closed" NCP, with 2 com- 
plete turns) and (3 from 0° to 360° (one period). A contour plot 
of our results obtained for h = 192 bp is shown in Fig. 0a). 
According to our definition, Lk(a, (3) also depends on the 
mean repeat length n trough the term P°/2jt = (« — 146) /7io 
{modulo 2k). Anyway, changing h simply shifts the levels of 
the contour lines without affecting their shapes. The last two 
terms of Eq. (|5} arise from the bending of the fiber. Taking 
into account that the mean radius of curvature of the fiber axis 
is equal to the fiber persistence length, which is ~ 30 nm fioll . 
it is possible to show that their net contribution per nucleo- 
some never exceeds 0.01. We can therefore conclude that a 
Lfe-conserving path in the (a, p) plane practically never devi- 
ates from a given contour line. 

We finally come to the issue of whether a biologically rel- 
evant Lfc-conserving decondensation path in the (a, p) plane 
can actually exist. We know that such a path must remain 
close to a contour line. Moreover, it should be capable of 
transforming a highly compact structure into a decondensed 
one. A fiber compaction map is displayed in color scale in 
Fig. [3j a). An exaustive scan of this map shows that there is 
only one small region, around point 1, which provides a de- 
gree of compaction high enough to match the maximal densi- 
ties observed in vitro, i.e. ~ 6NCP/10 nm (red in the color 
scale). Hence, there can be only one Lfc-conserving decon- 
densation path, namely the contour line going through point 1 
and leading to point 4. It is clear from Fig. [2 a) that varying 
only a (leading to point 2) or only (3 (leading to point 3) is 
forbidden at constant Lk. 

Moreover, energy considerations have to be taken into ac- 
count. As shown by Langowski in a recent simulation 12011 . 
the internucleosomal interactions play a minor role in the con- 
trol of the fiber compaction while the elastic energy (bending 
and twisting) of the linkers plays the major part. It comes out 
that the linkers should remain relaxed all over the deconden- 
sation path. First, in order to have relaxed linkers at point 1, 
P must be equal to p e<? (a°,n): this leads to select a special 
value for the mean repeat length h, equal to 192 bp (modulo 
10.6 bp). Secondly, during the decondensation (path 1^1), a 



varies from ~ 60° to ~ 90°. This variation corresponds to 
the wrapping of ~ 10 bp of the linker onto the protein spool. 
As the DNA helical repeat is lower in the NCP (h = 10.2 
bp/turn) than in the linkers (ho = 10.6 bp/turn), this wrap- 
ping is accompanied by a change of the equilibrium value 
$ eq (a,n) from 120°_to 140°. Remarkably, along the decon- 
densation path 1-4, p varies as well from 120° to 140°. Hence 
P = $ eq {a,h) all over the path 1^1, provided that h = \92bp, 
which is the repeat length measured in HeLa cells. 

Our results provide a new way of understanding the quanti- 
zation of the repeat lengths 12111 . The average repeat length 
has to be selected, e.g. thanks to nucleosome positioning 
sequences, in order to fix the mean nucleosome orientation 
p e<? (a,n). This particular orientation not only provides the 
maximal fiber compaction with relaxed linkers 12 ill , but also 
enables the global conservation of the linking number of the 
fiber during its winding/unwinding process driven by the nu- 
cleosome internal dynamics itself. 
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